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Abstract 

We propose to apply the two-particle irreducible (2PI) formalism to the problem of thermalization in 
heavy-ion collisions in the Color Glass Condensate (CGC) picture. We consider the 2PI effective action 
to three loops and derive a set of coupled equations for the classical Yang-Mills field and the quantum 
fluctuations in the boost invariant coordinate system. The initial condition and the relation to previous 
works are also discussed. 



1 Introduction 

Ultra-relativistic heavy-ion collision experiments at BNL-RHIC and at the CERN-LHC 
offer a unique opportunity to give us a glimpse of QCD matter under extreme conditions and 
the possible formation of the quark-gluon plasma (QGP) [IJ. This matter is intrinsically non- 
stationary, and understanding its expansion ('the little bang') bears a close parallel with tracing 
the history of the early universe at the dawn of the QCD epoch. Experimental observables 
related to collective flow and particle correlations put strong constraints on the late stage of the 
evolution which appears to be well described by hydrodynamics [2J. On the other hand, much 
less is understood about the early stage of the evolution, say, within 1 fm/c after the collision 
when the equation of state is not yet established, and one thus has to deal with the specific 
dynamics of the microscopic degrees of freedom in QCD. 

The color glass condensate (CGC) Q is arguably the most solid framework to date to de- 
scribe the very early stage of the nucleus collision when the coherent, strong color fields of the 
incoming nuclei are instantly liberated. What happens slightly after is vastly more complex. To 
first approximation, the fields of the CGC continue to evolve according to the classical Yang- 
Mills equation in the forward light-cone While this picture is useful for computing certain 
observables, it alone has little to do with the problem of thermalization. The latter inevitably 
requires the consideration of quantum fluctuations which, fueled by the decay of the classical 
field, grow explosively and drastically change the fate of the evolution. In the context of CGC, 
the importance of the fluctuations was first realized in [5], though they were introduced only 
heuristically, as random variables. A more proper treatment of the fluctuations as dynamical 
quantum fields has been implemented in II6I7I8II . The current state of the art [I7I8I9I is that one 
has a well-defined resummation scheme (originally devised for a scalar theory in which 
the classical solution is dressed up by quantum fluctuations to all orders. This will be reviewed 
in the discussion section. 
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Meanwhile, over the past decades significant progress has been made in the first-principle 
calculations of quantum field theories out of equilibrium, which goes under the name of the 
two-particle irreducible (2PI) formalism, or more generally, the n-particle irreducible (nPT) 
formalism [1 IJ. The 2PI formalism was initially developed for zero and finite temperature field 
theories, and is based on the Cornwall-Jackiw-Tomboulis (CJT) effective action r[0, G] [fT2ll 
which ab initio treats (the Green's function of) the quantum fluctuations G on equal footing 
with the classical field 0. When applied to the real-time evolution of nonequilibrium systems, 
the 2PI formalism can be naturally viewed as the field-theoretic generalization [il3J of the clas- 
sic, self-consistent ('^-derivable') method in quantum statistical mechanics II14I15II . A variety 
of nonequilibrium processes have been studied in this framework mostly in the context of scalar 
field theories HI 611711 8I19I20I2 1 1221231 . Li particular, it has been demonstrated that, starting 
from an arbitrary initial condition far from equilibrium, the 2PI dynamics drives the system 
towards the quantum equilibrium characterized by the Bose-Einstein distribution. On the other 
hand, there are not many applications to gauge theories due partly to subtleties in formal theory 
regarding gauge invariance and renormalization [I24I25I26I27II . So far, practical simulations of 
non-Abelian gauge theories far from equilibrium have been limited to the so-called classical 
statistical approximation [128129 ^ (see, also, f30l) which has the merit of allowing one to include 
2PI diagrams to all orders in the loop expansion at weak coupling, but systematically neglects 
at each order certain quantum contributions which are necessary to achieve the quantum equili- 
bration. The very same approximation is actually involved in the CGC-based approach [7 8] as 
we shall see in the discussion section. 

In this paper, we derive a set of coupled equations describing the evolution of the classi- 
cal color field a la CGC and its quantum fluctuations from the 2PI action to three loops in the 
boost invariant coordinate system. By truncating the expansion to fixed order, we do not include 
higher loop 2PI diagrams as in the classical statistical method, but we do not employ the clas- 
sical approximation, either. Our work is basically a marriage of the CGC approach and the 2PI 
formalism. While the latter has been developed, in part at least, with the motivation of studying 
the problem of thermalization in heavy-ion collisions, there does not seem to be a previous 
application of the formalism to the realistic setup of the collisions. We attempt to fill this gap. 
[See [|3T|| for an earlier work in the flat coordinates in the absence of the classical field.] 

The paper is organized as follows. In Section 2, we describe the initial stage of heavy- 
ion collisions as seen from the CGC viewpoint. In Section 3, we give a brief review of the 2PI 
formalism. We then derive in Section 4 the nonequilibrium evolution equations in the coordinate 
space from the 2PI effective action. The results are very complicated already at two loops, so in 
Section 5 we make a simplifying assumption and discuss the equations in the momentum space 
along with the initial condition. We then conclude in Section 6 by comparing our results with 
the previous CGC approach [17181 from a diagrammatic point of view. 



2 Gluodynamics in the A" = Q gauge 

The matter created in the central region of ultra-relativistic heavy-ion collisions at RHIC 
and at the LHC is almost baryon-free and nearly boost invariant. To first approximation, one 
may describe such a system as purely gluonic matter with strict boost invariance. A convenient 
choice of the coordinates is then the 'r-r^' coordinates defined by 
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r] = tanh ^ 
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x± = {x , X ) . 



(1) 



The 'proper time' r plays the role of ordinary time t, and r] is the rapidity. The components of 
the momentum in these coordinates are 



r 



Prj = X^p^ - tp^ 



P± = (P1,P2) • 



(2) 



Boost invariance means that observables are independent of 77, and the conjugate (dimension- 
less) momentum]?^ is conserved. An immediate consequence of this is that, around midrapidity 
where x^ ^ and t ^ r, the longitudinal momentum decreases as ^3 ^ prj/r. 

We shall be interested in the r-evolution of the gluonic matter in the forward light-cone 
r > 0. A popular choice of gauge is 



A'' = Ar = -{x-A+ + X+A-) = 

T 



(3) 



which is the analog of the temporal axial gauge A^ = in the ordinary coordinates. The Yang- 
Mills action in this gauge is 



Sym = J TdTdri(fx±_ 
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(4) 



where sums over the transverse directions i, j = 1,2 and the color indices a, 6, .. = 1, .., — 1 
are understood. We have lowered the Lorentz indices on fields, which we shall do throughout 
this paper. Introducing the notation = [r], x^) and the spatial metric 7q,/3 = diag (r^, 1, 1), 
the action can be compactly written as 



YM 



dTdri(fx^^ 



(5) 



where 7 = is the determinant of 7^^. 

The classical equation of motion following from the action ([5]) is 

= -^d. [V^l'-'d^A,) + r'l'^D.Fs, = , 

^ 6Aa ^/l ^ ' 

or in components. 



(6) 



-drA^ + ^D.Fir, = 



T \T / T 

--dr (rdrA,) + {\Dr,Fr,, + DjFji 



0. 



This should be supplemented with the Gauss's law constraint 

D^E"" = DiE' + D.^E'^ = , 
where the conjugate momenta E°' = yf^'-^'^^^drAp are 



(7) 



(8) 
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E' = TdrAi , E" = -drA„ . (9) 

T 

The energy density at the classical level is given by 

= ^{drA,r + ^{d^A,f + + ^4 . (10) 

However E = J Tdr](Px±£ is not constant in r because the metric 'ja/s depends on r. The 
continuity equation reads 

dr (^j Td7]d\^S^ = -T^ J dr]d^X±T'"^ , (11) 

where the 'pressure' in the //-direction is 

r'T^' = -^2(9rA,Y + ^{drA,Y + ^,F,,F,, - 1^ . (12) 

In the framework of the color glass condensate (CGC), the solution Aa of ^ represents the 
strong color field liberated from the colliding nuclei. At early stages, it is parametrically of order 
Aa ~ Qs/g where Qs is the so-called saturation momentum. Immediately after the collision at 
r = 0, it can be written as JH 

A = Al + A'i, A, = , (13) 
rdrA, = , = -gUcAl'A'i' , (14) 

T ' 

where A^ and A^ are the color fields of the projectile (nucleus 1) and the target (nucleus 2), 
respectively, before the collision. They are 'pure gauge' in the transverse plane diAj — djAi — 
ig[Ai,Aj] = and are related to the color charge density of valence partons p via diA]'"^ = 
p^'^. In the McLerran-Venugopalan model [32|, one calculates observables (such as 

the gluon multiplicity and energy density) using the solution A to the classical equation of 
motion which in turn is a functional of the initial fields A^''^. One then averages over A^'"^, or 
equivalently, p^'^ 

{O) = J Vp'Vp' W[p']W[p']0[A[p]] , (15) 
where the Gaussian weight functional 

W[p] ~ exp J d'^x±Pa{x±)pa{x±)^ , (16) 

with p ~ 0{Qs) takes into account the randomness of the color charge distribution in the trans- 
verse plane of a heavy nucleus. 
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Expanding the action around the classical solution Aa + a^, one finds the linearized 

equation of motion for the quantum fluctuations 



'1 

= 0, (17) 

where = 'y^'^Do.Dp and Da = da — igAa is the covariant derivative constructed from the 
classical solution. [Note that [Da, -D/j] = —igFap.] The interaction terms in the action read 

Sint = J dTdri(fx±^/j (C^ + £4) , (18) 
with the three- and four-point vertices 

A + A = -gfabc7''^l^\Daai,ra':^a's - ^^fabjab'c'r^l^'a'aalay; ■ (19) 
For later use, we introduce a shorthand notation 

Cbb',cc' = fabcfab'c' + fabc'fab'c , (20) 

which is symmetric under b b', c c' and (66') ^ (cc'). With this definition, (fT9l) can be 
written as 



-^a6',cc'7"Va''a^'<a^ (21) 



3 The two-particle irreducible (2PI) formalism 



The equations for the classical field Q and the fluctuations (fTTl ) by themselves do not contain 
the essential dynamics that drives the system to quantum equilibration. In order to go beyond, 
one has to consider the nonlinear interaction of the fluctuations shown in (fT9l) and its backre- 
action to the classical field. As explained in the Introduction, a powerful framework to discuss 
these issues is the two-particle irreducible (2PI) formalism which we shall review in this sec- 
tion. 

Real-time phenomena are best formulated in the Keldysh formalism which involves the dou- 
bling of the field degrees of freedom living on the closed time path (CTP) — a union of the for- 
ward [to, oo] and backward [oo, tq] time branches, [tq ~ is the initial time.] We define the 
propagator of the fluctuation 

G%{x,y) = iJc{al{x)4{y)}), (22) 
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where the symbol Tc denotes path-ordering along the CTP and the brackets (...) denote aver- 
aging over the density matrix initially prepared at r = tq. The effective action of the classical 
field A and the quantum fluctuation G is given by 

T[A, G] = Sym[A] + ^trlnG-i + '-ti G,'[A]G + T,[A, G] , (23) 

where r2 is the sum of 2PI diagrams constructed from the exact propagators (|22|) and the bare 
interaction vertices (fT9l ). Gq is the propagator in the presence of the classical background A and 
satisfies (c.f. (flTI) ) 

= s-s--J-S^^r - T')Sir] - r]')S\x^ - x'^) , (24) 
a/7 



where the contour-delta function is defined by 

M.-.') = K7"''; -'-',1^;-°^]' (25) 

[-6{t-t), t,t e [oo,ro]. 
The equations of motion for A and G are derived from the stationary conditions 
6T 

and 

To lowest order, (|26|) coincides with the Yang-Mills equation Q. But it receives quantum cor- 
rections (backreaction) from the fluctuation G as we shall soon see. Likewise, G deviates from 
Go due to the self-interaction of the fluctuations. (l27l) may be written in a more familiar form 
(suppressing indices) 

G = Go + GoUG , (28) 

where n(x,a;') = 2i6T2/SG{x, x') is the self-energy which has both the local (tadpole) part 
proportional to 5^{x — x') and the nonlocal part. Since the local part receives contributions only 
at one-loop and can be easily separated out, in the following we shall use the notation n(x, x') 
only for the nonlocal part. 

One of the merits of the 2PI formalism is the absence of the so-called secular terms IfTTTl . The 
secular terms appear ubiquitously in perturbative approaches to evolution equations in time, and 
have to be resummed by some means. Since the equations (|26l) and (l27l) in the 2PI formalism 
are those for the exact Green's function which resums infinitely many IPI diagrams, the secular 
terms do not arise at any stage. The price to pay is that one has to solve nonlinear equations. 
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(a) 



(b) 




(c) 



(d) 



(e) 



Fig. 1. 2PI diagrams up to three loops in the presence of the classical field. Crosses indicate the insertion 
of the covariant derivative. 

This is however straightforward to do numerically as an initial value problem, thanks to the 
manifest causal structure of the equations. 

In this paper, we assume weak coupling g <^ 1 and consider r2 up to three loops. The 
relevant diagrams are depicted in Fig. 1. In these diagrams, crosses represent the insertion of 
the covariant derivative D[A] = d — igA where A is the classical field. At early times when 
A ~ 0{l/g), the two terms in the covariant derivative are of the same order in the coupling. In 
this regime, the evolution is driven by Diagram[ll5 which describes the decay (decoherence) of 
the classical field, the scattering of quantum fluctuations off the classical field and the 1^2 
processes of the fluctuations. [The last process was previously considered in [31 J in the flat 
metric] With increasing r, the magnitude of A decreases, and so do the diagrams with insertions 
of A's. When A ~ 0{g^), they become of the same order as the three-loop diagrams. 

Comments are in order concerning gauge invariance. We first assure the reader that the resid- 
ual gauge symmetry in the A'^ = gauge will be respected in the evolution equations to be de- 
rived below, in the sense that they transform covariantly under residual (r-independent) gauge 
transformations. However, there are potential problems regarding more general (r-dependent) 
gauge transformations: By truncating the loop expansion at a fixed order in the coupling, one 
cannot guarantee the full gauge invariance of the 2PI effective actionQ This is a well-known 
problem of the 2PI approach and the general solution is not known. Still one can rely on the 
idea of "controlled gauge invariance" [124125 1 which states that the gauge dependence appears at 
higher order than the truncation order for the propagator, and at twice the order of truncation for 
physical quantities (pressure, etc.). Moreover, in QED there exist generalized Ward identities 
for the correlation functions which save renormalizability [I26I27II . It would be very important 



We note, however, that the problem at hand (gluodynamics in the t — r] coordinates) has never been 
discussed outside the A'^ = gauge to our knowledge. 
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to study these issues in the present context. We leave such a study for future work. 

Before leaving this section we make a technical observation which will greatly simplify the 
subsequent analysis. The equations shown in the previous section are cluttered with factors of 
the metric tensor 70/3. One can eliminate them by introducing a dimensionful coordinate C = tt], 
and accordingly, the following rescaling 



Aj^ = tA(^ , ttr, = ra^ , 9^ = . (29) 
The equation of motion (|7]) then becomes 



7(M,)/r 



r 



1 5Sym (pa ' 1 



- d'^A, + -drA, - D^F^, - DjFj, . (30) 



SA, V " r 

Note that one should not do the rescaling (|29l) before functionally differentiating the kinetic 
term of the action because it interferes with the r-derivative. Let us label the new set of spatial 
coordinates = {(, x_l), di = {d(, d±) by capital letters /, J, ... By slight abuse of notation, 
we denote the two equations in (l30l) as 



- d^^Ai + -drAj - 5,^4 - DjFjj , (31) 



keeping the above caveat in mind. Eq. (|3TI) suggests that in the coordinate system {(, x±) the 
equations (l26l) and (ITT]) look as if we were working in the Cartesian coordinate system. Indeed, 
the interaction vertices (|32|) take the same form as in flat space 



Sint = J^dT£x(^-gfaUdiajTa'}a'j - ^—Cw,cc'A^ia'ja]a'j - |-a6',cc'a/ «j«7«j) >(32) 

where d^x = d(d'^x_L and the subscript C means that the r-integral is performed along the CTP. 
Similarly, the equation for the Green's function (fTTI) becomes 



-1 ab 



dl + -dr 



T 



Sij - {D^5i.j + DiDj - 2DjDj) 
5jK5'''Scir-T')S%xr-x'j), 



(33) 



where D'^ = DjDj and we have also rescaled the Green's functions G^^, = tt'Gcc, Gni = tG 



'CO 



Girj 



r'Gii^. Since (|32|) does not involve r-derivatives, no subtlety arises when taking the func- 
tional derivative of diagrams constructed from these vertices. Thus the use of the coordinates 
x^ leads to the simplest derivation and representation of the evolution equation without any 
complications from the metric tensor. The only cautionary remark is that one has to remember 
which value of r is used in the rescaling C, = tt] in each of the covariant derivatives that 
appear in the equations. 
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4 Nonequilibrium evolution equations 



In this section we explicitly write down the evolution equations (|26l ) and dTT] ) using the two- 
loop effective action. The contributions from the three-loop action are complicated and thus 
relegated to the Appendix. Actually, at the level of the derivation of the evolution equations we 
do not have to assume boost invariance, and the results in this section are valid for generic back- 
grounds Aj depending on all coordinates {t,'1],x±). Assumptions on the (r^, xx)-dependence 
come later, when we specify the initial conditions and make approximations to solve these 
equations. 

4. 1 Evolution equations in the coordinate space 

Let us now evaluate the functional derivative in (|26l) . Consider first the term 

^ (^trGo^G) ~ i^(tr {D%j + DjDj - 2DjD,rG'j^^ , (34) 

where the symbol ~ means we only keep terms which contain A. Explicitly, 

itr {D'rc'fj f dx (di^iAlGfj) + A\dj^GJj + gf^^^AlAiGj;) _ , (35) 

hi {DjDjrGf, ~ ^ X {diiA'jGfj) + A'.djGfj + gUA'jAjGfj)^^^ . (36) 
Differentiating, we obtain 



- -tr {DTG'!j = yUd^Gfjix, y) - d^Gfj{y, x 



V2 ' ' ■'V 2 V ^ -^-^^ '''' ' 'Jy=x 

'\' — {fbacfcde + fbdcfcae)AjGjj{x,x) 

= pabcD':j{Grj{x, y) + GJAy, x)),=. , (37) 
where = d^5^'^ + gfbdeAj{x) and similarly, 

^ ^hr {DjDjrG'f,) = p,,,{d^jGfj{x, y) - d^jGfj{y, x)) 



SAf V2 



^2 ^ ^ 

'^'^{fbacfcdeAjGjj + fbdcfcaeAjG^jj) 

p,,,D':j{Gyjix, y) + Gfjiy, x))y=, . (38) 



Noting that y) = Gfj{y, x), we find 

6 (% 



V2 



(^tr GtG) = gfabc [olWfAx, y) + D':j(GJ,{x, y) - 2Gr,(x, y 



(39) 



y=x 
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(a) 



(b) 



(c) 





(d) 



(e) 



Fig. 2. The self energy diagrams up to two loops (from the three-loop effective action). 

relevant diagram is Fig. [ij) with the interaction 



Next we calculate 5V2/5A. To two loops, the 
vertices shown in (l32l) 



contractions 



= -^Ca,,c, I dy VLn.LMNG'k{^. y)GTMi^, y)G%{x, y) , 
where we defined the (dressed) three-gluon vertex 

LMN — —ig^^LMC^flm'nD^ + flmri'D^"^) 



(40) 



(41) 



In (|401) and in the following, it is understood that each covariance derivative in V acts only on 
one of the G"s which has the same color index. Combining Eqs. (|3TI) . (|39| ) and (l40l) . one finds 
the equation of motion for the background in the presence of quantum fluctuations. 



D'^jGyjix, y) + D':j(Gy,{x, y) - 2G- (x, y) 

I dyVl^^,^,^Gf,{x,y)GTM{x,y)G%{x,y) 
J c 



9 fabc 

f 

2 



9 ^ 



(42) 



At the three-loop order, there are additional contributions, (lA.ll) and (IA.2I) . to the right hand 
side of this equation. 

The equation for the Green's function (l28l) reads 
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1 ah 



fr + - ^^ic) ^iJ - (DHij + DjDj - 2DjDj) 

= -i6''''6jK6c{x-y), (43) 

where we have separated out the local part of the self energy coming from Fig.[T^. To one-loop 
(two-loop in r2), the nonlocal part comes from a single diagram Fig. [2^ and is given by 

n/j(a;, y) = - '^almJLAlGLL'i^^ y)G^[lii{x, y)VbVm',JL'M' ■ (^4) 

To two loops (r2 to three loops), there are five more self energy diagrams Fig. [2b-2e which 
should be added to (|44|) . These are evaluated in (|A.3I) - (IA.6I) . 



It is important to notice that the above equations (1421) and (1431) transform covariantly under 
the residual (r-independent) gauge transformation of the background 

Ai UAiU^ + -UdiU^ , (45) 
9 

and the concomitant transformation of the fluctuations 

a? ^ {UaiU^Y = U^^a] , (46) 
GfAx^y) ^ UfGi^\x,y){Ulf'' = UfufGif{x,y) , 

so that the equations take the same form after a t/-transformation[l] 
This can be checked by using the following identifies 

f T rbb' T TCc' T raa' e f T T<m' T rbb' t tccI f 

Jabc>J <J — U Ja'b'c' , Jabc<J U U — Ja'b'd , 

Cab^cdU'^'U^^'U''' = U''^'Ca'y,c'd' , (47) 

which hold for a S\]{Nc) matrix U'^^ = (t/^)^" = (U^^)'"^ in the adjoint representation. In fact, 
the two terms of the interaction Lagrangian (fT9l) are already separately invariant under the trans- 



^ The term —A(^/t'^ on the left hand side of (|42] ) may seem problematic at first. However, in the trans- 
formation 

Ac UAcU^ + -UdcU^ = UAcU^ + --UdM^ , 

g gr ' 

the second term depends on r, and the inhomogeneous terms cancel among the first three terms of (l42l) 
so that 



as it should. 
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formations (|45] ) and (|46l) which are a decomposition of A + a — [/(A + a)W + ^UdW . Hence 
the derivatives of diagrams constructed from these vertices are expected to behave properly 
under the (residual) gauge transformation. 



4.2 The statistical and spectral functions 

In practice, the equations (l42l) and (l43l) are awkward to handle because of the complicated 
time integral defined along the CTP. A standard trick to disentangle this is to rewrite the equa- 
tions in terms of the so-called statistical J' and spectral p parts of the Green's functiord] 

G{x,x') = {Tc{a{x)a{x')}) = ^({a,a'}) + ^{edr - r') - Qcir' - r))(Ka']) 

= \^Qc[r - r') - Qc{t' - r)) p . (48) 

Roughly, the spectral function tells which states exist, while the statistical function tells how 
many particles are there in each state. The spectral function is related to the retarded Green's 
function Gr as 

Gr{x, x') = -i9{r - t')p{x, x') . (49) 
In the same way, we define the statistical and spectral parts of the nonlocal self-energy 

n = - '-{edr - r') - 0c{r' - r)) Up . (50) 

From the definition ([48]) in terms of the (anti-)commutator of fields, it follows that J-'{x, x') and 
p(x, x') are real. That njr(a;, x') and np(x, x') are also real is not immediately obvious, but can 
be understood from the identities 

^f = GrJ^Gj^ , 

^p = KG-r'-Ga^). (51) 

where Ga{x, x') = i9{T' — t)p{x, x') is the advanced propagator, and the fact that Gr and Ga 
are purely imaginary in the coordinate space. 

It is then a simple exercise to show that (|42|) can be rewritten as 



^ In the literature the statistical function is often denoted as F. Here we use the calligraphic letter in 
order to avoid confusion with the field strength tensor Fjj. For the spectral function we use the same 
letter p as the source color charge (c.f., ([T5])). but the distinction should be obvious from the context. 
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- gUc (D':jJ^Zi^, y) + D':j{rfj{x, y) - 2J--(x, y) 

/ y=x 

d'y yy,^ U (x, y)FJ^,{x, y)J'%{x, y) 



,^9 ^ 

'T—^yyab.cd 



TO 



(52) 



where d^y = T'dT'dr]'(fy± and the r-integral is now defined normally. [For clarity, we omitted 
indices in the last three terms in the large brackets. They are the same as in the first term. We 
shall use similar abbreviations below.] In the d^y integral one may prefer to use the retarded 
Green's function instead of p 



d^yp{.x,y) 



d'^yiGR{x,y) 



Similarly, (1431) decomposes into the real and imaginary parts 



(53) 



r 



+ -d^ - -6j^ 6jj - {D'5ij - DjDj - 2tgFjj] 



ab 



^TA^.y) 



- r d'z U,{x, z)fjT'fj,{z, y) + r d'z n^(x, z)fjp%{z, y) , (54) 



fr + - ^5i^) Su - {DHjj - DjDj - 2igFu) p'^x, y) 

+9^ (Cad,beJ^Ijix,x) + \Cab4e^MMi.X.x)^Ij] P^jA^.y) 

\ab J)C 



ab 



d''zU,{x,z)fjp%{z,y), (55) 



where 



^^F\X-,y)ij ~ 2 almJLM \^ LL'-' MM' ^PlL'PmM' ] ^ bl'm',JL'M' 



bl'm',JL'M' ' 

xy 



(56) 



and 

np(a^, y)u = '^almJLM {J^LL'PmK-P + PlL'-^MAp) ^bl'm'JL'M' ■ (5^) 



It is straightforward to include the three-loop contributions to r2 in the above set of equations. 
However, the result is rather cumbersome and we do not show them in this paper. 

In principle, given the initial condition one can solve these equations numerically and study 
the degree of thermalization and isotropization by looking at various observables. [Observables 
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should be averaged using the weight function as in (fT5l).1 In practice, however, this might be 
difficult because of the complicated structure of the equations already at two loops. In order to 
render them amenable to numerical simulations, in the next section we consider the special case 
where the classical field A is homogeneous. 

5 The case of the homogeneous background 

5. 1 Evolution equations in the momentum space 

To solve (|52l ). (|54l) and (|55] ) in full generality is a daunting task as it is computationally too 
expensive. The main source of difficulty is the spatial inhomogeneity of the background field 
Aj. [If we take boost invariance for granted, inhomogeneity here means the dependence on 
x±.] As a matter of fact, in the 2PI literature the problem of solving nonequilibrium evolution 
equations in the presence of inhomogeneous backgrounds has not been tackled numerically 
even in scalar field theories. We thus assume from now on that the background is homogeneous, 
leaving the inhomogeneous case for future work. From the viewpoint of the application to CGC, 
admittedly such an assumption is unrealistic because it is in general incompatible with the initial 
condition ([H) Ai = A} + Af being dependent on x±, dtA]' = p^'^ 7^ 0. Nevertheless, since the 
equations are quite nontrivial even after this simplification, one may expect that they still capture 
the dominant features of thermalization and collective behaviors averaged over distances larger 
than 1/Qs , the typical scale of variation of A. 

If the background is homogeneous, we can define the spatial Fourier transform of Zu = 
{J", p, Ujr, Up}ij as 

Zjj{r,T',7^-r]',x^-x'^) = J -0-^e^P^(v-v')+^P.i^.-^'JzJJiT,T\p), (58) 

where pa = iPr],P±) and d^p = dprj(fp±. With this definition J-'{p) and p{p) are dimensionless 
while njr(p) and Iip{p) have dimension 4. Note that we cannot perform the Fourier transform 
with respect to the rescaled momentum o C. Nevertheless, the notation pj = (p^, p±) is still 
convenient. The covariant derivative acting on the left of the Green's function becomes, in 
the momentum space, 

(9c - igA^)I{r, r', v-r]') = -(9, - tgA^)Z{T, r', r/ - r]') 

r 

-{ip^ - igA^)Z{T, t\p^) = {ipc^ - igA(^)Z{T, t\p^) . (59) 
r 

In the following, we use pa in loop integrals d^p and in the argument of the Green's function, 
whereas we use pi in the covariant derivative keeping in mind which value of r is used to rescale 

PC =Vrilr. 

In scalar theories with a homogeneous background, Z(r, r', p) is real due to parity. However, 
in gauge theories this is no longer the case. Since implicitly depends on A, parity symmetry 
and the condition that Z is real in the coordinate space 

Z{p, A) = Zi-p, -A) = A) , (60) 
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only require that the real part of Z(p) is even in p while the imaginary part is odd in p. Thus, 
for instance the first term in the second line of (|52l) becomes 



idf6'^^ + gf,,eA'})TTjix,y) 
d?p 



y=x 



d?p 
(27r)3 



{iPi5''' + gUA%T)):Ffj{T,T,p) 



(2 

Similarly, in the cubic terms of (l52l) 

1"^ I, I f d^p d^q 



I {-Pi'^^^fAr, T,p) + gUA'Ar) TleTfjij, T,p)) . (61) 



X 



1 



pUr, r\p)7TM{r. r\ q)F%{r, r\ -p - q) + {FpF) + {FFp) - -{ppp) 



(62) 



with the dressed gluon vertex in the momentum space 



V,. 



lmn,LMN 



{p,q,^) 



-igyLM^iiPN - qN)flmn + gClm,npAlf{T 

+SLN{i{-'2pM - qM)flmn + gCln,rapA\j{T 

+SMN[i{PL + '2qL)flmn + gCmn,lpAl{T) 



(63) 



the terms in V which are linear in momenta do not vanish even though the background A does 
not carry spatial momentum. 

The equations for the fluctuation (|54l) . (|55l) become manifestly complex-valued 



1 



T T 



+ 



-{tPLS"" + gfaedADizpLS'" + g U,Ai)^',^^{T, t',p) 
{ipid'-' + gfaeciA]){tpj5''' + ghfbAi) + 2g^UMhefA]AQ F%{t, t\p) 



(27r)3 V^"'^'^"'^"^^' ^' ^ 2'^-'^'deJ'MM{^^ ) ^' .P) 



= - rT"dT"U,iT,T",p)fjJ^Ur",r',p)+ r T"dr"UAr,r",p)fjp%{T",T',p), 

J TO J TO 

where the self-energy in the momentum space is 

n^(r, r",p)fj = 1 1 ^T/_n,/M7v(p, q, r) (j-]S^',(r, r" , -q)J'TN,{r, r",p + q) 

-^Pmm'(^, -(1)p7n'{^^ r",p + q)^Vbm'n',jM'N'{-p, -q, t") , 
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+pMip{r, r", -g) JT^,(r, r",p + q)]v,^,n',jM'N'{-V. -q, r") . (64) 



5.2 The initial condition 

Having discussed the evolution equations in the momentum space, we now specify the initial 
conditions at some initial time r = tq ~ 0. The initial conditions for the classical field are 
already stated in (fT4)) . The second equation in ([14)) implies ~ C(r^) as r — 0, from which 
the initial condition for = A^/r follows 



= 0, drA^ = -^U^Al'A^' . (65) 



On the other hand, the initial conditions for the fluctuations J-'{tq, tq, p) turn out to be nontrivial 
and have been a subject of recent debate. Commonly in the 2PI literature, J-'(ro, tq, p) is chosen 
to be an arbitrary function far from the equilibrium distribution (like the Gaussian or 'tsunami' 
distribution). However, in our problem the initial condition is not arbitrary, but fixed by the 
underlying QCD Lagrangian. 

Let us first consider the free theory without the classical background. The initial distribution 
of the fluctuation ('Wigner function') derived in |l6l is equivalent, in our terminology, to the 
initial condition for the statistical part of the Green's function. The results of Q translate into 



•^"(ro, ro,p) = 5., i/., [ + ^ ^.^^..r^ , (66) 



drdr'^F^, (r, r ,p) = S^b^ 5,, - , , ~ ^abS'^,-^ , (67) 

where the last expressions are valid when | ^ tq |p_L | . As for the initial condition for J^^^, one 
should note that the longitudinal field a,, = ra^ and their correlation functions are constrained 
by the Gauss's law 

-d„dra„ = dnC^ = —Did = —rDidrai , (68) 
r 

as obtained by linearizing ([8]) around the classical background E"^ + e\ —^E^ + e^. 

Using T^fj = tt'J^c^c^ ~ 0(r^r'^), one finds, in the absence of the background, 

-^ct(^0'^o,p)~0, (69) 
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4 p2 



2 2 

,„,^=£^_^5,,-^, (70) 



where we assumed 7^ 0, see below. 

The presence of the background Ai will modify the above formulae. Roughly, pi — )• Pi — gAj. 
More precisely, as discussed in [8J, one has to diagonalize the matrix M in the quadratic fomo 



-{DLDLT'Si, + {DiD.Y' + 2igF^' a) , (71) 



for each realization of A, from which one finds ~ A^^^/^. Clearly, this procedure will affect 
only the subleading terms in Tq so that for our purpose they may be neglected by taking tq 
arbitrarily small (except in (|70l)). For instance, if we discretize the spatial coordinates (r/, x^) 
by putting them in a periodic box of size (L.^, L^) = (2Nnari, 2N±aj_) where a,^ and aj_ are the 
lattice spacings, the condition ^ To|pj_|, ro|Aj| reduces to 



2vr|r2^| 



Taking ^(Aj ~ Qs < l/aj_ and 2Nj^a^ ~ 5, the typical size of the (approximately) boost 
invariant region in heavy-ion collisions, we see that the above condition can be well satisfied 
even for tt,^ = 1 if we use a very anisotropic lattice tq ^ a^. Such an anisotropic lattice is 
commonly used in the 2PI literature. As for the zero mode = = 0, we can neglect it 
altogether at initial time since is then independent of rj and can therefore be absorbed into 
the initial background A^ which itself is a random variable in CGC. We thus conclude that the 
leading (tq — > 0) terms of (|66l ) and (|67] ) may be used as the initial condition of our equations in 
the presence of the background, whereas (|70l ) is modified by the Gauss's law (1681 ) to 

= - ^9facaA!i){-5''p, - ighe,A^^ ^p\ + tIp\ - 

^ {5'"'p^ - igfacdA-r) {6'% + ighedAI) . (73) 



4|P.| 

For r > To, the evolution equations automatically keep terms of order r^pf and T'^{gAiY. 

Finally we derive the initial conditions for the spectral part of the Green's function. The 
canonical commutation relation (in the momentum space) 

[ai, rdraj] = iSij , [a^, -9^a^] = i , (74) 

T 



^ If the background A depends on x±, this involves the diagonalization of a large matrix with the di- 
mension proportional to the spatial volume. 



17 



immediately gives 



Being derived solely from the canonical commutation relation, Eq. ([751) actually holds for all 
values of r, and moreover, even in the presence of the background field. The former property 
will serve as a good check of the numerical simulation. Eq. (fTSl) together with the equation of 
motion for p^^ (as r — )■ 0) and the antisymmetric property Pr;»?(T, r') = —p-rirjir', r) uniquely fix 
the small-r behavior 

Pvvir, r',p) ^ \{r^ - r'2) , (r, r' ^ 0) . (76) 
One thus finds the initial condition for p^^ 

PCc(r, r',p)^ ^4^^ , (r, r' ^ 0) , 9,pg(r, r,) = -5^' . {11) 
2rr To 



6 Discussions 

An efficient and illuminating way to describe the nature of the equations (|52l) . (|54l) and (1551) 
would be to compare with the previous CGC-based works. The approach pursued in [7 8] is 
an extension to QCD of a technique developed for a scalar theory ifTOll . One first solves the 
Yang-Mills equation @ perturbatively 



Ai = Af^ + AY' + AY' + 



(1) 



l(2) 



aP ~ 0{g 



n-l\ 



(78) 



where Af^ is itself an exact solution to the Yang-Mills equation and A^' satisfies the linearized 
equation 



(1) 



ij 



DfDf 



2^gFf:^ 



A 



(i)b 



0. 



(79) 



The initial condition is taken to be A{tq) = yl(°)(ro) + A^'^\to). Since A^^^ + A^^^ is not an 
exact solution, one finds the (linearized) equation of motion for A^"^^ 



dl 



T 



IJ 



Df^Df 



-I ab 



A 



{2)h 



(80) 



Note that the right hand side is identical to the second line of (I52l) if one replaces y) — )■ 
J-'^'^\x,y) = A^^\x) A^^\y) . Let us write (f80l) schematically using the retarded Green's func- 
tion G^-* oc p^^^ in the background of A'^^^ 



9P 



(0) £)(0)yl(l)^(l) 



(81) 
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(a) 



(b) 



^ >^ ^ • 

(c) 



t 



Fig. 3. Tree diagrams obtained by the perturbative expansion dTSl ). Crosses and blobs denote the insertions 
of D^^^ and A^^\ respectively. Solid lines connecting vertices are the retarded propagators oc p^^^ 
in the background of A^^\ 

and represent this as in Fig. |3^. Repeating this procedure, one finds two contributions 
which are depicted by the two diagrams in Fig. [3b, respectively, and 



(82) 



(83) 



corresponding to the diagrams in Fig. [3]:. Continuing in this way, one generates an infinite 
number of tree diagrams made up of D^^\ A^^^ and the retarded propagator. One then makes 
'contractions' of pairs of A^^^ in all possible ways (by using an appropriate Gaussian weight 
functional of A^^\tq)) and identify each pair with J'^'^^ As a result of this, one gets loop dia- 
grams as depicted in Fig.lH [Note that A^'^^ vanishes after this operation since it has an odd num- 
ber of A'^^^'s.] Each building block of these diagrams has a counterpart in (|52l) . (|54|) and (l55l) : 
Diagram 4b corresponds to the terms ~ pJ-'J-' in (|52|) . The tadpole on the retarded propagator in 
Diagram 4c represents the modification of the p-function which is taken into account by the J-'p 
terms in (l55l) . The structure p(o)j^(o)p{o)p{o) in Diagram 4d corresponds to the IIpp ~ {^p)p 
terms in (|55]), p(°) J'^o) J'^o) in Diagram 4e corresponds to the J" J" terms in ([541), p(°) p(°) J'^o) J'^") 
of Diagram 4f corresponds to the TipF ~ (pj-") J-" term in (|54l ), and p*^°) J^^^^ J^*^°)p*^°) in Diagram 
4g corresponds to the IIj-p ~ {J^^)p term in (|54l ). 

Higher order diagrams for A*^^'^'- which are solely made up of the above building blocks 
are automatically resummed by the nonlinear equations (I52l)-(l55l). However, diagrams which 
contain more complicated irreducible structures can only be resummed by evaluating r2 to three 
loops (as in (|A.3I) - (|A.6I) ) and beyond. These higher order contributions are initially suppressed 



19 



(a) 



(b) 



(c) 



X X — t-. X X 



I 

X -X 



_ ' 

(d) (e) (f) 



(g) 



Fig. 4. Loop diagrams obtained from Fig. |3] after performing all possible contractions of all the blobs. 
Dashed lines denote T^^^ . 

by powers of the coupling A^") ~ but at late times they may become important since they 
are multiplied by factors of J-' which may grow at low momentum due to some sort of instability 
(see, e.g., [l33l and references therein). 

On the other hand, already at two loops there are terms in (|52l) - (|55l) which have no coun- 
terpart in the above diagrammatic expansion. These are the ppp term in (|52|) and the pp term 
in Iljr. [There are no loops formed by two retarded propagators.] Neglecting these terms is 
precisely the content of the classical statistical approximation ifTTI . Thus one can view the 
CGC-based approach Il7l8ll as a generalization of the classical statistical approach II28I29II to 
the situation where there is a nonzero background: In the former, one characterizes the large 
occupation number of gluons by the (inhomogeneous) background A'^°)[p] ~ 1/g associated 
with the source charge p, while the initial distribution of the fluctuations J-'(ro, tq) at tq ~ 
is uniquely fixed by the QCD Lagrangian expanded to quadratic order. In the classical statisti- 
cal approach, one usually sets A^^^ = assuming there is no source (or by gauge invariance), 
and characterizes the large occupation number by an arbitrary, but physically motivated func- 
tion J^(r°,r°) = (A^^^tq) A^^^tq)) according to which the initial configurations A^^^tq) are 
generated. This provides an effective description of the evolution after the classical field has 
decayed substantially, meaning that tq > l/Qs- 

As is well known, the classical statistical approximation is valid only for low momentum 
modes where the occupation number is large. In a free massless scalar theory in flat space. 



where n{p) is the Bose-Einstein distribution at the temperature T. The condition | ^ |p| is 
valid for momenta much smaller than the temperature p <^T so that n{p) ^ 1, and even in this 




(84) 
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regime the approximation is valid only in some averaged sense due to the oscillation in t — t'. 
Because of this, the classical approach cannot reproduce the Bose-Einstein distribution, or its 
Boltzmann tail. Instead, one finds [,34.9J 



n{p) 



T 



1 



2 



(85) 



P 



as the equilibrium distribution as determined from the vanishing condition of the collision in- 
tegral HpJ^ — Ujr p (c.f., (l54l)). The two terms in (f85l) are precisely the first two terms in the 
expansion of the Bose-Einstein distribution at p ^ T. 

In contrast, in the 2PI formalism one is guaranteed to obtain the Bose-Einstein distribution 
as the equilibrium distribution because by including the pp terms, or rather, by not neglecting 
any term one can satisfy the detailed balance within each of the diagrams contributing to 11. 
This in particular means that the formalism encompasses a relatively large (ideally, an infinitely 
large) region in momentum reaching out to the Boltzmann (exponential) distribution which 
is an unmistakable feature of heavy-ion collisions. In principle, the present approach has the 
potential to describe the evolution of the system from right after the collision all the way to the 
late quantum regime in a single framework. 

In conclusion, we have set up the foundation to incorporate the physics of the CGC in the 
2PI formalism. As already remarked in Section 3, there still remain important formal problems 
concerning gauge invariance and renormalizability II24I25I26I27II which have to be addressed in 
order to place this approach on firmer ground. And of course, eventually one would like to solve 
the evolution equations numerically and apply the results to the phenomenology in heavy-ion 
collisions. At least in the homogeneous case, the causal structure of the equation allows for a 
straightforward numerical implementation (provided some renormalization scheme) as an initial 
value problem. Yet, manipulations of large matrices such as J^ij{t,t' ,pa) may pose serious 
challenges for limited computer resources. Work in this direction is under way [[351 . 
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A The 2PI action to three-loops 

In this appendix we consider r2 to three loops and calculate its functional derivatives. There 
are three diagrams as depicted in Fig.[lt-[It- The last two diagrams give the following contri- 
butions to be added to the right hand side of ([421) . 
Diagram [T]l: 
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^p(ld) _^2 r ^ , ^ 

J^(^ = ^^^b,im J dydzdwVlp,^ pp,jiG%{x,y)G^lp{x,z)G'J^^{x,w) 

^^i[_R'(^) '^)G^P'L'{^i y)^'Q'M'h^i y)^qq'r',QQ'R'^b'l'm',I'L'M' ■ (A.l) 

Diagram [l^: 



X 



2/)Gi^p(a;, + 2G%,{x, y)G%{x, z) 

X + 2Gi«Q(x,z)G^r;(^,y)) ■ (A.2) 

The two-loop self-energy diagrams coming from Fig[l]:-[It are shown in Fig.[2b-2e. They 
give corrections to (l44l) and are evaluated as 
Diagram 2b: 

{Jl^^Y^j{x, y) = -^Caa' ,U'Cw ,mra'G^J^j^{x, ?/) ( j G^m + 2G"^ G^j ) . (A.3) 



xy 



Diagram 2c: 



{Ii''')fj{x, y) = \j dzdw \^:im,ILM^;p'r,PP'RGUx, ^)G7q(^, ^) 

^^RR'i^-i '^)GpiLi{z, y)GQ,]^j,{w, y)'^^g'r' ,QQ'R''^bl'm' ,.JL'M' ■ (^-4) 



Diagram 2d: 



(n^^)/j(a;, y) — —^Cppi^qqi j dzV^ii, jj;^j;^,V,^^i jj^j^pG lp{x, z) 

X (g{U^, z)Gi^t'{^, y) + 2G%{x, z)G%{z, y)) y) . (A.5) 



Diagram 2e: 

■ 2 

(n^'')7j(a:, = ^Cbb/^„„/ j dz \^llm,ILM^pqr,PQRGLp{^^ ^) 

+(a/ ^ 6J) . (A.6) 
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